Entanglement Entropies in Conformal Systems with Boundaries 
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We study the entanglement entropies in one-dimensional open critical systems, whose effective 
description is a conformal field theory, both for the ground state and the excited states associated 
to primary fields. The analytical results are checked numerically finding excellent agreement for the 
c = 1/2 minimal model and the c = 1 compactified massless free boson. 
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In the last years it has been recognized that quan- 
tum entanglement is the key concept to understand the 
intricate structure of many body wave functions in Con- 
densed Matter Physics. This recognition led to the in- 
vention of novel numerical methods, such as MERA and 
PEPS, which generalize the DMRG method (see [1] for a 
review). The success of these numerical techniques lies in 
the correct implementation of the area law that governs 
the entanglement entropy (EE) of the low energy states 
of local Hamiltonians [2] • In one spatial dimension, the 
entropic area law implies that the EE is bounded: this is 
the typical situation for systems with a gap in the spec- 
trum. However, if the system is gapless and described by 
a conformal field theory (CFT), the EE grows with the 
log of the subsystem size which amounts to a violation 
of the entropic area law. It is well known that the dom- 
inant finite-size correction of the ground state energy of 
critical systems is related with the central charge c [3]. 
A similar correspondence also exists between the central 
charge and the Renyi entanglement entropies (REE's) of 
the ground state (GS) @HS]- In the same vein, we would 
expect that the scaling dimensions of the primary oper- 
ators are related with the Renyi entropy of the excited 
states (ES's), since the finite mass gaps of critical lattice 
Hamiltonians are also related with them. Indeed, this 
latter correspondence was proven recently for conformal 
systems with periodic boundary conditions [3 [8]. The 
main result of these works is that the n th Renyi entropy 
is given in terms of the 2ro-point correlator of the cor- 
responding primary field placed at special positions that 
depends on the subsystem size. 

In this letter we shall generalize this result to open lat- 
tice Hamiltonians described by boundary conformal field 
theory (BCFT) [§]. BCFT has a wide range of applica- 
tions to impurity problems [10] . string theory [11] . etc, 
which justifies the generalization pursued here. More- 
over, BCFT has a very rich mathematical structure which 
is worth to analyze from the entanglement point of view. 

We shall start with some basic definitions. The REE's 



are defined as follows: 
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where n is a positive real, pa is the reduced density ma- 
trix of a subsystem A and the trace is over the A's Hilbert 
space. If pa is constructed from the GS of a CFT, one 
has EJ: 
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where I is the size of A, L is the total-system size and r\ = 
1,2 for periodic/free boundary conditions (PBC/FBC); 
c is the central charge and is a constant that depends 
on the BC's [T2]. For open systems, I is measured from 
the left edge. Equation ^ may have significant correc- 
tions which carry useful information about the under- 
lying CFT [TSTfTB] . Different kind of corrections arise 
when, instead of considering the GS, one looks at the 
ES's [16l [17], or when the system satisfies general 
BC's preserving its conformal invariance (FBC are just 
one of them) [TS] [TS] . In particular, the corrections in 
this last case have, up to now, not have been derived an- 
alytically, and their knowledge would be a great advance, 
both from a practical and a conceptual point of view. In 
the following, we present the general CFT framework al- 
lowing the analytical computation of such corrections. 
These results are then verified with DMRG calculations 
in two examples: the c = 1/2 minimal CFT and the c = 1 
massless compactified free boson. 

Let us consider a ID conformal system, defined on a 
strip of width L. Cardy [9] showed the existence of BC's, 
denoted as {a}, that preserve the conformal invariance, 
and such that the partition function takes the form 
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where q = e % P/ L (being /3 the inverse temperature), h 
are the conformal dimensions of the primary fields, and 
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Xh are the associated characters. The integers Af£p are 
the fusion coefficients of the theory [T2l [20] . The deriva- 
tion of the REE's of ES's, originating from the primary 
operator T(w) for open conformal systems follows closely 
the one of periodic boundary conditions given in [7J . The 
REE's are given by 
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where pa,t and Pa,gs are the reduced density matrices 
obtained from the ES produced by the primary field T 
and from the GS, respectively. The ratio can be 

computed using path integral methods [7J. We explain 
below the main steps of the derivation. 

Let us first split the strip into the subsystems A = 
[0, l],B = [I, L]. This strip is parameterized by the com- 
plex coordinate w — a + ir, a g [0, L], t G (— oo, oo). We 
make the conformal transformation 
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which maps the strip into the unit disc D = {£, |£| < 1}. 
The intervals A and B are mapped into the segments 
(—1,0) and (0,1) respectively. The boundaries of the 
strip, a = 0, L are mapped into the boundary of the disc 
(|C| = 1), and the infinite past (w^ — —ioo) and infinite 
future (w^q — ioo) are mapped into 
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oo ^oo 
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Next, we make n copies of the unit disc and sew them 
along the cut (—1, 0), obtaining the Riemann surface lZ n , 
which can also be mapped into the unit disc by the con- 
formal transformation 
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The points give rise to 2n points on 7Z n 
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where the primary fields T and its conjugate are in- 
serted. Repeating the same steps as in [7J, we arrive at 
an expression for Eq. Q 
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where we have assumed that T is a chiral primary field 
with conformal dimension h. A similar formula holds 



for non chiral fields. The correlators in Eq.(|9| are com- 
puted on the unit disk for fields inserted at its boundary. 
These primary fields change the boundary conditions a 
and b on each edge of the strip: T changes them from 
a to b, while from b to a. Hence the numerator of 
Eq. (|9| is proportional to the partition function of a disc 
with 2n segments where the boundary conditions a and 
b alternate [12] . Equation ([9| constitutes the main result 
of the work, which we shall verify below for two models. 

The c = 1/2 minimal CFT. This CFT contains three 
primary fields: the identity I (hj = 0), a Majorana 
fermion \ (h x — 1/2) and a spin field a (h a = 1/16), 
whose fusion rules are 
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This CFT describes the long distance properties of the 
critical Ising model in a transverse field whose lattice 
Hamiltonian is 



H, 
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where a x ' y ' z are the Pauli matrices. The correspondence 
between the conformal BC's {a} and the lattice BC's is 

the following: (|), corresponds to fix af L to +1 (-1), 
while corresponds to free BC's [HI HE]- For simplicity, 



we shall denote 0, |, ^ by 



,F. 



We begin by considering the REE's of the GS with 
FFBC's (the two F's refer to free BC's on both edges). 
The spin Hamiltonian ( |11[ ) can be mapped into a free 
fermion Hamiltonian |2l|. which we use to compute, by 
using the method of [22] . the von Neumann entropy, i.e., 
the n = 1 REE, for chains with lengths 60-180 in multi- 
ples of 20 (this will be the sizes considered in the rest of 
the work). Using Eq. and finite-size scaling (FSS) 
techniques, we obtain the asymptotic values c = 0.499 



and 



„<7=2 _ 



0.241. The value of c agrees to high precision 



(7) with the central charge of the critical Ising model. More- 
over, the value of c\~ 2 is very close to c\~ x /2, that can 
be obtained with PBC. The relation c\~ 2 — c r [~ 1 /2 = ln g 
O [TB] , is satisfied in this case, because the boundary en- 
tropy (BE) In g vanishes for FF BC's [5S]. Moreover, as 
a different reliability check, we verified that the REE's 
in this case are exactly one half of the ones of an XX 
spin- 1/2 chain with FF BC's (see below), according to 
reference pM] . 

We next study the +- l-BC case (that is equivalent to 



BC). The fusion rule 7\/J 



d 1 ' 



lead us to consider 



only the case where h = 0, for which no corrections arise 
(apart, as we shall see, from constant BE contributions). 
We compute the REE's for the Hamiltonian (11) using 



the DMRG method [53] with up to 800 states per block, 
and 3 sweeps, which yields a truncation error of 10~ 12 or 
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FIG. 1: REE's in the ID critical Ising model. Black to 
turquoise (dotted lines): L — 60 to 180 numerical data; pur- 
ple (solid line) + arrow: CFT predictions (the BE's are added 
to the CFT formulas when necessary). 



The Figs. [ljc) ,(d) show the convergence of the numerical 
results to the CFT prediction obtained with (12 1, for 
n = 1, 2. 

We then consider the +FBC case. The fusion rule 
Aff a = 1 implies that the spectrum contains just the 
conformal tower of the a field, so the F function must be 

(n) 

F„ . The nontrivial fusion rule of the field a yields dif- 
ferent chiral correlators, which were computed for general 
n in [28] . In our case though, only a combination of them 
yields the appropriate F functions, which for n = 2, 3 are 
given by 
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The Figs. [jje), (f) display the DMRG results and the 
CFT predictions: the agreement, for L — > oo is excellent, 
up to the constant term — | In 2. 

Finally, we consider the H — BC case (equivalent to the 
— hBC case). The spectrum contains just the conformal 
tower of the x operator, and therefore the REE's of the 
GS shall be the same of the ones of the first ES in the 
FFBC case (see Figs. [TJg) , (h)). To conclude, we have 
shown that the CFT predictions for the Ising model with 
all possible conformal BC's agree with the numerical re- 
sults, up to the constant contribution due to the BE. 

The c = 1 compactified free boson. We shall next con- 
sider a massless free boson with a compactification radius 
R = 1. This CFT is rational, meaning that the chiral 
symmetry is enhanced in such a way that the number 
of primary fields is finite and have conformal dimensions 
h = 0,1/8,1/2,1/8 [H]. The corresponding conformal 
characters are given by setting A = 0, 1, 2, 3 mod 4 
(h x = A 2 /8) in 



Kx(q) 



(14) 



less. The Figs. [TJa), (b) display the results for n = 2, 3 
REE's: the data progressively flattens to the theoreti- 
cal value -\ In 2 PH |2B]. The convergence to the CFT 
predictions is of order 10 -3 or less, as confirmed by a 
FSS analysis. This behavior can be ascribed to the pres- 
ence of a slowly F-depending finite-size correction [24], 
previously observed, for the Ising model, in [5]. 

We now consider the .FFBC case. The fusion rule 
N£ a = 1 implies that the first ES is generated by the 
primary field x- The associated F function is given by 

where 



F- 



(n) 



is a massless free boson with 



x 

c = 1 [Sj. Quite interestingly, F^ has a general expres- 
sion valid for any value of n > 27J 
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where 77(g) is the Dedekind function. The conformal BC's 
are Dirichlet (D) and Neumann (N), and the correspond- 
ing boundary entropies are and -|ln2 [TU]. 

The lattice realization of this CFT is given by a spin- 
1/2 XX chain, with boundary couplings [29] 

L-l 

#B = -E( CT M+i + CT Kfi) + 
1 

_ 2 



(15) 



where the DBC, on a given edge, is realized by setting 
all the boundary couplings to zero, while the iVBC is 
realized, say on the left side, by choosing a z = and 
a + = a_ = 2 (moreover, in the sermonic picture of the 
F>£>BC case, one has to work in the half-filled sector). 
The operator content of the various models described by 
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the Hamiltonian (15) can be expressed in terms of the 



partition functions of a free boson with different BC's 
[TO] 129] . The DD and NN partition functions can be 



written in terms of the characters ( 14 1 



Z DD (q) = K (q)+K 2 (q), 
Znn ((?) = K (q). 



(16) 
(17) 



The absence of corrections for the GS in the DDBC case, 
with the exception of the usual oscillating ones [15] , has 
already been observed |30| ) and we confirm them numer- 
ically. We expect the same feature for the GS in the 
NNBC case, which we analyze with DMRG for system 
size L = 100, keeping up to 1100 states, using 3 sweeps 
and achieving a truncation error of 10 -10 . The results 
are shown in Figs. [2](a), (b), for n = 2,3 : up to os- 
cillating corrections, typical of c = 1 systems, and the 
BE — | In 2, as expected from Eq. (|9j) we do not see any 
non-constant correction. 

We then consider the first ES with DDBC which, ac- 
cording to Eq. (16), is generated by a vertex opera- 



tor with conformal dimension 1/2. The F^ 1 ' functions 
of vertex operators T are equal to 1, for all n, so the 
REE's receive no corrections T]- We verify this result 
by identifying this excitation in the fermionic version of 



the Hamiltonian ( 15 ) with the addition (or subtraction) 
of a fcrmion at half filling. The Figs. [2jc), (d) show 
the n = 1, 2 REE's for these states, which confirmed the 
absence of corrections. 

For DD and NNBC, there is another ES, with con- 
formal weight h = 1, associated to the primary field 
idcf>. For DDBC it corresponds to the lowest particle- 
hole state created from the half-filled ground state [7J [5] ; 
in the NNBC case, there is no particle number conser- 
vation, so it is simply the first ES. In the former case, we 
use the method of [22], and in the latter the multi-target 
DMRG [3T], to compute the relative REE's, shown in 
Figs^e), (f), (g), (h). Up to oscillations (and the BE 
in the NN case), we find excellent agreement with Eq. 
@. 

Finally, we study the NDBC case, whose partition 
function [TU] cannot be written in terms of the characters 
(14). However, it can be shown that 



z N D(q) = xi/w{q)[xo{q) + Xi/2(q)], 



(18) 



where Xo,i/2,i/i6 are the characters of the primary fields 
of the Ising model. Equation ( 18 1 implies that the op- 



erator content corresponding to NDBC is given by the 
tensor product of two Ising models. In particular, the 
GS is associated with the operator 1(g) a, resulting in the 

(n) 

correction F„ , and the first ES with a ® \i resulting in 
the correction Fa Fx (plus a BE contribution — \ In 2): 
we show in Fig. [3] the results for the n = 2, 3, obtained 
with multi-target DMRG, finding, even in this CclSG, a 
remarkable agreement between CFT and numerics. 
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FIG. 2: REE's in the model |15| (I). Black dotted lines: L = 
100 numerical data; red solid line: CFT predictions (the BE's 
are added to the CFT formulas when necessary). 



We conclude that for the critical Ising and XX quan- 
tum spin chains under any conformal BC's, the REE's of 
the low-lying states are obtained, apart from the oscilla- 
tions and the BE's, from the correlators of the primary 
operators of the underlying CFT: the finite-size behav- 
ior of the REE's of quantum chains with open bound- 
aries identify the primary operators in the CFT. Since 
the REE is simple to evaluate with numerical methods 
like DMRG, MERA and PEPS, these results, beside in- 
teresting, should have a large appliability. 
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FIG. 3: REE's in the model ( 15 1 (II). See the caption of figure 
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